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Abstract 

We investigate mean-field effects in two- component trapped Fermi gases in the superfluid phase, 

lol we 



in the vicinity of s-wave Feshbach resonances. Within the resonance superfluidity approach 
calculate the ground state energy and the heat capacity as function of temperature. Heat capacity 
is analyzed for different trap aspect ratios. We find that trap anisotropy is an important factor 
in determining both the value of heat capacity near the transition temperature and the transition 
temperature itself. 
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At present several groups 



have produced Fermi gases of atoms at temperatures 
where the superfluid phase is expected. The possibihty of tuning interatomic interactions 
near Feshbach resonances may provide a chance to watch macroscopic phenomena for both 
normal and superfluid phases and the phase transition from one to another. The presently 
reachable temperatures limit the existence of the superfluid phase to a strongly- interacting 
gas. For now it is still an open Question how to observe the superfluid transition, although 
numerous proposals exist (and references therein). 

Here we examine the thermodynamics of a two-component Fermi gas in a strongly- in- 
teracting regime where we can expect superfluidity. Although any attractive interaction can 
theoretically support Cooper pairing at zero temperature, it seems that only a strongly- at- 
tractive interaction may produce observable effects in current experiments. For a strongly- 
interacting gas the 'small parameter' kpa > 1 and so the ordinary mean-field approach 



breaks down. In this case the resonance superfluidity [10|, [111, [1^ approach suits the prob- 
lem best. On the experimental side the strongly- interacting regime is reached by exploiting 
magnetic-field Feshbach resonances For our numerical simulation we took the 

parameters of the two- component Fermi gas of '^^K near Feshbach resonance. These com- 
ponents are in|9/2 — 9/2> and |9/2 — 5/2 > states. In 1] it was shown that there is an 
s- wave Feshbach resonance at a magnetic field 224.21 ± 0.05G and near this resonance the 
scattering length was measured up to about ±2000ao. 

KokUelmaus et a. Q employed the resonance superfluidity approach in order to get the 
critical temperature and find possible signs of the phase transition. Ohashi and Griffin |12| 
independently used the same approach. For the superfluid phase we use the mean-field ap- 
proach of BCS-type theory within the local density approximation and use the Hamiltonian 



of the resonance superfluidity model ll| in order to take Feshbach resonance physics into 
account. In this case the energy density functional S depends on normal p(x, y), anoma- 
lous k(x, y) densities and a molecular field where x and y are coordinates of two 
atoms. The set of equations as well as their solution for energy- independent effective in- 
teraction(see for example) look like in ordinary Bogolyubov-de Gennes formulation but 
with the extra term of a molecular field: 

p(k, R) = n(k, K)u^{k, R) + (1 - n(k, R))f ^(k, R), 
nik, R) = u(k, R)v{k, R)(l - 2n{k, R)) 
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E(k,R) = ^/i(k,R)2 + A(R)2, 



\v^{k,R)J 2' ^(k,R)' 

h{k, R) = ^ + Vmf{R) + VtrapiR) " A (l) 



W) 



A(R) = -Vtg I ^^/€(k,R) - Q(j)(R) 



where k is a wavevektor, R = is the coordinate of the center of the mass for two 
interacting atoms, n(k, R) = (exp(-E'(k, R)/A;bT) + 1)^^ is the Fermi-Dirac distribution, 
Vmf{R) is a mean field potential and A(R) is the energy gap. V^g = Anh^ahg/m, ai,g is the 



background scattering length, g = ^VtgABAjj is the coupling strength and, u = {B — Bq)Ah 
is the magnetic field detuning, AB is the field width of the resonance, A/i is the rnagnetic 



moment difference between two hyper- fine levels of the two-component Fermi gas 
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—-^d'R{ep{k, R) + HpP(k, R)p(R) 



More complete theoretical analysis using this approach can be found in 12|. Functions 
f{k,R) are the Wigner transforms of corresponding functions f{x,y). So after solving 
we have both a normal p(k, R) and an anomalous ^(k, R) distribution function as well as 
the molecular field 0(R) and we find the ground state energy as: 

_fk_ 
(2^ 

+p(k, R)2Vtrap{n) - A(R)/t*(k, R)} (2) 
+ J d'RVZtim'i^) 

It is a well- known phenomenon that near the phase transition the energy-temperature 
curve should have a distinct change. Moreover we can investigate the dependence of specific 
heat capacity and ground state energy on interaction strength and trap aspect ratio. As we 
will see the geometry of the trap strongly influences the thermodynamics of the superfluid 
gas. 

We consider an anisotropic trap with trap aspect ratio A = ujz/u!± of the transverse 



and axia 



0.01) 



m 



frequencies. Because the ratio A is rather small in the current experiments 



calculations for isotropic or almost isotropic traps possess only a methodological 



interest. However, it is in nearly isotropic traps that the most dramatic observables are 



expected; see below. As the interaction near a Feshbach resonance strongly depends on 
detuning we consider all the characteristics of this system as a function of detuning v and 
temperature T . For present calculations we chose lo^ = AOOHz and 5 x 10^ particles. As the 
detuning is defined by the magnitude of a magnetic field we will report detuning in units of 
Gauss, which is more convenient for a possible comparison with experiment. 

Within this approach we have calculated the ground state energy and heat capacity for 
a variety of interactions and temperatures. The typical dependence of ground state energy 
and gap on detuning is shown in Fig.( HE))- In the range of detuning shown the interaction 
is always attractive so a superfluid phase can exist. But at a detuning larger than ^ IG 
the pairing energy seems rather small compare to the kinetic energy as the energy- detuning 
dependence is very weak. In the vicinity of a resonance the mean-field Vmf contribution to 
the ground state energy is negligible as the background part of the scattering length is small 
compared to the resonance part. With decreasing detuning the pairing energy becomes 
a significant part of the total energy and we can see considerable lowering of the total 
energy beginning from some detuning that depends on temperature too. Unfortunately it is 
impossible to estimate the kinetic and the pairing energies separately because the integral 
of the second and forth terms in (0) over momentum are individually divergent Q]. 

Fig.^thus shows that detunings smaller than ^ IG seem required to observe superfluidity. 
In Fig.Q]the solid curves are calculations for A = 0.5. The dotted curves are for A = 0.05; 0.01 
and for T = O.lTp. So for more anisotropic traps the pairing energy is smaller. It should 
be possible to measure the energy, as in In the released energy was measured for 
Bose gas above and below the transition temperature. The uncertainty is ~ 10 — 30% for 
temperatures below T^. The desired detuning should therefore be small enough to generate 
at least a 10% change in energy compared to the normal phase energy. Thus the detuning 
should be smaller than 0.5G for T /Tp = 0.1, A = 0.5. For more realistic trap aspect ratios 
(0.01) it should be smaller than 0.3G At these values of detuning the pairing energy is not a 
small part of the total ground state energy which means that the energy gap is rather large 
too (Fig. El). 

From experimental point of view it seems that the closer to the resonance we are the more 
chances to observe superfluidity we have. But it should be mentioned that near a resonance 
there is so called the BCS-BEC crossover regime 3]- In theory of the superconductivity at 
this regime there is a class of 'exotic', high-Tc superconductors. In the case of the crossover 
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regime (/i/A < 1) we need a more sophisticated theoretical approach than the ordinary 
mean- field theory. Such a treatment is outside of scope of the present article, but it may 
be important in order to get the right value of the critical temperature [2^ for this regime. 
If /i/A is sufficiently greater than 1 we can hope we have the ordinary BCS regime. For 
smaller trap ratios the condition (/i = A) happens at a slightly smaller detunings. 

In this article we do not analyze how much larger /i/A should be for the theory to 
hold, and leave such an analysis for the near future. So for further investigations we have 
chosen detunings v = 0.3 — 0.5G where the energy gap is still significant. The energy- 
temperature dependence (Fig. IH]) is a more appropriate characteristic for the comparison 
and describing of experimental data. Again from an experimental uncertainty point of 
view it seems that temperatures around O.lTp or smaller are more appropriate to observe 
superfluidity. For chosen parameters of interest we have the condition A ^ hu and according 
to the pairing classification of this is an inter-shell pairing regime for which the local 
density approximation is the appropriate tool for describing the pairing. 

The important characteristic of the pairing field is its distribution in momentum 
n(k, R) (Fig.lH). In the same figure we demonstrate this distribution for u = 0.4; IG (dashed 
and dotted curves) for T = O.lTp. Far from resonance {u = IG, dashed line) this function 
has just a narrow peak near Fermi momentum kp- With decreasing detuning more and 
more atoms below kp are involved in the pairing which emerges not only as a Fermi surface 
effect. As it was shown in the product kp^ (^- the coherence length) is the appropriate 
variable for 'exotic' superconductors. Moreover the value ^ can be considered as the size 
of the Cooper pair. Within our approach the coherence length is dependent on location 
^^(R) = / (irK(r,R)r^/ / drK{r,K) ^ {kp(K)/rmiA(R))'^. We have calculated this for the 
center of the trap and found that ^ is about 0.33-1.65 in trap units for u = 0.3 — 0.5G and 
A = 0.01 But the interparticle distance is changing much slower and is around 0.35 for the 
given detunings. For our reasonable detunings and A = 0.01 the parameter kpC, is 2.26-8 



but the 'boundary' which distinguishes the high-Tc and conventional superconductors [2l| is 
kpC^ ^ 271. It was suggested [3] that the 'exotic' superconductors are intermediate between 
BCS- type supeconductors and BEC. Holland el al IQj predicted the same phenomenon 
for the two-component degenerate Fermi gas. So according to this classification we have 
conventional superfluidity for u ^ 0.5G and larger. Moreover at this field kpa < 1 and the 
ordinary mean-field theory can be used. 
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The coherence length can give us an insight into whether it is possible to detect a signal of 
the Cooper pair breaking, similar to what was done for bound molecular states at negative 
magnetic- field detunings {2^. If it is possible, this can be considered as a sign of superflu- 
idity. From this point of view it seems that it is not desirable to be very close to resonance 
because in this case the many-body wave function of the Cooper pair at distances of the 
coherent length order between atoms will be very close to the two- body wave function for 
quasi- bound states near a resonance. At large distances between atoms the wave function 
of the Cooper pair will be very different from the relative wave function of two scattering 
atoms but it is nor clear if the RF spectroscopy can work at such distances. Our approach 
enables us to find the Cooper pair wave function and in the future we will analyze this aspect 
in more detail. 

The phase transition to superfluidity stipulates a considerable increase in heat capacity 
near the critical temperature. This characteristic may therefore serve as an observable sign of 
superfluidity. The main point of the present article is an investigation of a heat capacity for 
Fermi gas for the conditions described above. In the authors calculated and analyzed this 
characteristic for off- resonance interacting ^Li within ordinary HFB theory and the effect is 
too small to be detected in current experiments. But the resonance superfluidity approach 
gives us a chance to do this analysis for strongly- interacting systems. The phenomenon of 
superfluidity strongly depends on the trap geometry. For an s-wave interacting gas the R— 
dependence on the energy gap reflects only the trap geometry . In our calculations we 
keep uj± = 400Hz and vary Uz- It is clear that with decreasing A the energy gap will be 
smaller in the center of the trap than for an isotropic trap. So the spatial manifestation of 
the superfluid phase jll| depends on geometry too. 

We found (Fig. El EI) that the trap anisotropy considerably influences the heat capacity 
of a gas in the superfluid phase and can wash out the effect altogether if A is very small. It 
is known that at the critical temperature the heat capacity has a discontinuity. In the trap- 
confined gas the Tc as well as the gap is R- dependent and the higher the temperature is, the 
smaller is the fraction of the gas involved in superfluidity. The critical temperature is the 
point at which the energy gap disappears in the center of the trap. Above this temperature 
the gas is in the normal phase everywhere. We can see that for temperatures lower then 
Tc heat capacity of the superfluid phase can be up to ^ 1.5 times larger than for normal 
phase gas, and this value depends on the trap aspect ratio and detuning. Also the more 
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anisotropic a trap is the smaller the gap is and then smaller a 'bump' in heat capacity for 
superfluid phase. 

The heat capacity as well as the en ergy of ground state can be measured using a ballistic 
expansion as was done for a Bose gas [T^- But, as it was pointed out in an expanding 
strongly- interacting gas is likely in the hydrodynamic regime. We have estimated that for 
temperatures near the critical temperature and at small detunings (0.2-0.4G) the collision 
rate |2| can be larger than the trap frequency, which also depends on the trap aspect 
ratio. It means that the hydrodynamic regime is very probable during the expansion. We 
can choose the trap aspect ratio in order to get the collision rate a bit smaller than the 
transverse frequency. For u = 0.3G at the corresponding critical temperature the ratio A 
should be around 0.01 (Fig. E)). A similar estimation for u = OAG shows that A can be around 
0.05(Fig. ini). It is clear that for reasonable parameters near the resonance the cloud will 
not be expanded ballistically as a whole. However, a measurement of the total energy after 
turning-off the trap is still valid, since the released energy should be conserved. The trap 
energy is almost two times smaller than total energy so the energy- temperature dependence 
of the released energy will be very similar to that shown in Fig. IHl but approximately two 
times smaller. 

This work was supported by the NSF. The author appreciates J.L.Bohn and M.Holland 
for stimulating and useful discussions and C.Ticknor for providing with parameters of the 
resonance for "^^K. 
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FIG. 1: Ground state energy per particle as a function of the detuning for different temperatures. 
The temperature is in Tp units. The trap aspect ratio A = 0.5 for the solid curves. The dotted 
curves display different aspect ratios for T = O.lTp 
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FIG. 2: Energy gap in the center of the trap and chemical potential as a function of the detuning 
for different temperatures. A = 0.5. 
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FIG. 3: Ground state energy per particle as a function of a temperature for the magnetic field 
detuning v = 0.2 — 0.6G with step = O.IG The dotted line is the energy in the case when the 
interaction is described just by the non- resonant scattering length. 
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FIG. 4: Distribution function K;(k, 0) in the center of the trap versus momentum for the pairing 
field for different temperatures in the case v = 0.2G and A = 0.01. Also shown are the cases for 
u = 0.4(dotted line) and v = lG(dashed line) and for T = O.lTp. 
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FIG. 6: Heat capacity as a function of a temperature for detuning u = 0.4G and for different trap 
aspect ratios A. The dashed hne is for the normal phase gas. 
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